New method for the quantum ground states in one dimension 
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A simple, general and practically exact method is developed to calculate the ground states of 
ID macroscopic quantum systems with translational symmetry. Applied to the Hubbard model, a 
modest calculation reproduces the Bethe Ansatz results. 
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Since the very beginning of the quantum theory, to 
solve the Schrodinger equation for macroscopic quantum 
systems has been one of the main tasks of theoretical 
physics. It would not be an exaggeration to say that, 
due to lack of such methods, a considerable effort of the- 
oretical physicists has been devoted to the development 
of a variety of perturbative and approximate methods 
and numerical simulations. But a desire for powerful 
non-perturbative methods has grown stronger over the 
last couple of decades with the list of phenomena played 
by strongly correlated electrons getting longer, particu- 
larly since the discovery of high temperature supercon- 
ductivity in copper oxides [lj]. While we have seen a 
considerable progress in rigorous treatment of quantum 
ID and classical 2D systems over the last several decades 
0-0], these rigorous methods are not flexible enough to 
solve non-integrable models in one dimension, nor, most 
probably, generalizable to higher dimensions. On the 
other hand, the method of NRG (numerical renormal- 
ization group), particularly DMRG (density matrix RG) 
has seen a remarkable success first in quantum ID sys- 
tems fiol ] and then in finite Fermi systems, competing 
well with the conventional quantum chemistry calcula- 
More recently, the notion of entanglement 
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tions 

from quantum information theory 12j helped a further 
progress in NRG towards the finite temperature as well 
as dynamical quantities flijjl5j ]. 

In a recent article, we have developed a simple, general 
and practically exact method to calculate statistical me- 
chanical properties of macroscopic classical systems with 
translational symmetry up to three dimensions [lfjj]. We 
here extend this method to solve the Schrodinger equa- 
tion for ID quantum ground states with translational 
symmetry. As a benchmark model for this development, 
we consider the Hubbard model. Just like our recent 
work on the 3D Ising model, our method is purely alge- 
braic and other than seeking a convergence in entangle- 
ment space, it does not employ any other notions such 
as NRG, nor make any approximations. Our results for 
the ground state energy and the local magnetic moment 
in the ID Hubbard model agree with the known exact 
results by Bethe Ansatz 0, Qjjj ■ An important difference 
of the present method from the Bethe Ansatz, however, 
should be emphasized: the new method is not rigorous 



but mathematically much simpler, general and therefore 
readily applicable to any quantum spins, fermions and 
bosons. This is a reflection of the fact that our recent 
method for the Ising model is applicable to any classical 
statistical systems with translational symmetry. Yet an- 
other but probably the most significant remark here is 
that the success in ID Hubbard model should constitute 
an essential ingredient in the analysis of the 2D Hubbard 
model by the present method. Again, this is a reflection 
of the fact that our recent method for the 3D Ising model 
crucially relies on the successful analysis of the 2D Ising 
model, we called it the "Russian doll" structure, and the 
mathematical structure involving the D=2,3 Ising models 
and that for the D=l,2 Hubbard models are essentially 
identical. 

The Hubbard model is defined by the Hamiltonian, 



H = —t 
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where t is the tranfer integral, a measure of kinetic en- 
ergy, U is the onsite Coulomb potential and Cj CT , c\ a are 
the annihiration and creation operators for electrons at 
site i and spin a. We take t as the energy unit. To 
calculate the ground state of the Schrodinger equation 



(2) 



we follow the following steps. 

First, instead of ([2]), consider the eigenvalue problem 
for the density matrix 
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A well-known observation about ((3|) is that, starting with 
a trial wavcfunction ^ which has non-zero overlap with 
the ground state, only the ground state survives in the 
limit (3 —> oo. Monte Carlo and NRG simulations are 
based on this observation [HI [H. Here our idea goes 
opposite, j3 — » 0, and calculate the largest eigenvalue of 
the operator 1 — /3H and corresponding eigenstate. 

Second, we rewrite the Hamiltonian |lj as a sum of a 
local bond Hamiltonian, 



H=J2(H ij +H i + H j ) = 



bond 



bond 



bond 
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with 



H 



H 



a,<ij> 



h.c. 



U 
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where the onsite Coulomb term is split into two sites 
i and j, and the chemical potential /i is introduced to 
control the electron number per site. 

Third, we note a decomposition of the density matrix, 
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This is the simplest Suzuki- Trotter decomposition [19j, 
but it is good enough for (3 — > 0. In (J7]), following the pro- 
cedure familiar in quantum Monte Carlo, we have split 
the entire bonds into two groups: one connecting the sites 
(2i, 2i + 1), the even group, and the other (2i + 1, 2i + 2), 
the odd group. Now the local bond density matrix should 
be further decomposed as, 
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where and below the repeated indices imply a summa- 



tion, and Q a takes five operators, 1- 



(n 



if 



H,4.J, c,-^, Cif, c^, and C4 and Q a likewise operators at 
site j. Since the local pair density matrix contains 
even number of creation and annihiration operators, the 
matrix representation of the density matrix ([7]) can be 
written as a operator product of local matrices, 



(lk\ 
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%d \ij) ~ fa,ik ® 9a ,jl 
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where 



/1 = 5i 
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FIG. 1: Schematic figure of the transfer matrix eigenvalue 
equation, EqQ. 



etc., where four basis states at each site are ordered as 
|0), I |), II) and I tl)- Note that the -1 in the f 2 matrix 
is due to the fermion anticommutation algebra. Thus the 
matrix product representation of the even group bonds 
in the density matrix is, 



• fa ® 9a ® fp ® gp ® / 7 <S> g-y ■ 



(10) 



and the same expression for the odd group bonds with 
one lattice shifted from the even group case. Putting to- 
gether, we have the matrix representation of the density 
matrix (|7|) as, 



K 



g a ■ fp ® / 7 • gp ® g 1 ■ fs 1 
fe-gs®ge- fv®--- 

j5 1 
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where for notational simplicity, we have raised the in- 
dices 1,2 for the two Ts to their shoulders. Note also 
that F^g are 4x4 matrices for each pair of interaction in- 
dices (a, (3). Thus, r 1,2 are a set of 5 2 x4 2 numbers which 
will be denoted below like r^ 2 ,^, where (0,6) indicates 
(up, down) interaction channels, whereas (c, d) indicates 
(left, right) basis states. 

Fourth, we write the ground state wavefunction as, 



' ' ' (apm ® C| 7 a 2 ® C 7 «5a 3 ® C 



2 
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on the basis • • • |oi) ® |a 2 ) ® (03) <8> \<h) ® ' ' * where ai 
etc takes 4 states |0), t)i II) an d I tl)- One can derive 
the form (fl"2")) by a successive use of matrix algebra [lj| . 
Consider, for example, a wave function ^(ait^o^c^). Re- 
garding this as a matrix of the left index a± and the 
right index {a 2 a 3 a4}, SVD (singular value decomposi- 
tion) gives *(aia 2 a 3 a 4 ) = J2 a A a ia p a B{ a2a3a4 y a . The 
quantity B can in turn be regarded as a matrix of 
the left index {020;} and the right index {0304}, thus 
SVD gives B {a2a3a4}a = J2p c {a 2 a}p^pD{a, 3 a 4 }p- Like- 
wise, £>{a 3 a 4 }/3 = YlfEiasPh^Fa^. Putting together, 
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rewriting A aia as A a (a 1 ), C{a 2Q }/3 as C a p(a 2 ), -E{a 3 ^} 7 
as Epjiaz) and F ail as F 1 (aA), and appropriately ab- 
sorbing p a , Xp and A 7 into the matrices A, C, E and F, 
one gets ^(a^a^) = A a (a 1 )C a p(a2)Ep 7 (a 3 )F 1 (a 4 ). 
For our density matrix with a bipartite structure with 
translational symmetry, (jlip . one arrives at the claimed 
form. Again for notational simplicity, we have raised 
the indices 1,2 for the two £s to their right shoulders. In 
quantum information theory, these indices a, f3 and 7 are 
known as entanglement [12]. Considering only 1 for these 
indices is a simple mean- field- like approximation for 
Allowing larger values, one takes into account the effect of 
correlation with increasing precision. An important note 
here is that (fT2")l is not peculiar to the Hubbard model, 
but rather a general statement for macroscopic quantum 
ground states with translational symmetry. Putting the 
above arguments together, the eigenvalue problem ([3]) 
for (3 — > is then written schematically as in Fig. 1. The 
horizontal lines indicate 4 local basis states, whereas the 
vertical lines indicate 5 interaction channels connecting 
nearest neighbor sites for r * and entanglements for £ ' . 
To emphasize the close similarity to our recent analysis 
of the Ising model, let us call all 4 lines associated with 
r 1,2 as bonds. Note that Fig. 1 is a slight generalization 
of Fig. 1 in the 2D Ising model [16j ■ 

Fifth, we follow the procedure in our method for the 
Ising model, namely we handle the eigenvalue problem 
j3| as a variational problem. We thus maximize the 
quantity fj,Q = ^K^/^f^! by iteration starting with 
an input state for ^F First consider the numerator. 
A local ingredient of this quantity is, Au> etmm if = 

(L a T l gaa >(t>n> a >(n m p T2 gfP0>(n> m >p>- The real nonsym- 
metric matrix A can be written as A = RvL tr where the 
matrices L, R, and v are made up of left eigenvectors, 
right eigenvectors and eigenvalues of A and tr means the 
transpose. The eigenvectors are normalized as L tr -R = 1, 
and due to this property, the summation over the com- 
bined entanglement-bond indices We in the numerator 
can be done N — 1 times, N — > 00 in the end, and thus 
we only keep the largest eigenvalue vq and eigenvectors 
L and R . We have fAf = v$- x luj r AELq. The de- 
nominator is handled likewise. Let us denote the cor- 
responding largest eigenvalue and eigenvector as p§, Lo 
and Ro. Note that \iq now contains C 1 ' 2 in a quadratic 
form. Maximizing this quantity with respect to Q- and 
( 2 then leads to generalized eigenvalue problems: 

<S{L (mm'/)Ro [ll'e)T) gaa , T 2 gef)p , £ l/3 ^ n , Vfi , }Ci„ a 

= fi S{Lo(mm')R (ll')5 aa ,5 l 3i3>Cw(n'i'i3'}(Ln a (13) 

= ^{LoCmmORofMO^OQ'^/S'CmTiaCm'n'Q'}^ ( 14 ) 

where the symbol S means a matrix symmetrization and 
Mo = fk^o^/po' 1 - We solve (T3J) and (H]) for the next 
£ 1,2 and continue until convergence. 




-0_8 I 1 1 1 I 1 1 1 I 1 1 1 I 1 1 1 I 1 1 1 

0.0 0.2 0.4 0.6 0.8 1.0 

Electron number per site 



FIG. 2: The ground state energy per site at 17 = 8 vs. the 
electron number per site, n = 1 corresponding to half-filling. 
From the top, the entanglement n=l,4,6 (rectangles) and 8 
(star, red online). The thick solid line (blue online) is the 
Bethe Ansatz result [l7t ]. 

Finally after the convergence, we can calculate various 
ground state properties. In the present case, quantities 
of interest are, the average number of up spins (rif), that 
of down spins (n±), the double occupancy (nfU^), the 
kinetic energy per site — t J2cr{ c ia c j>? + h.c), and the lo- 
cal magnetic moment ((^f?) 2 ), where a is the Pauli spin 
matrix. In general, the expectation value for the two op- 
erators A and B sitting on the adjacent sites is calculated 
as. 

(AB) = ^AB*/** 

= {a2a 1 \AB\aW2)C0aA'P'a' 1 CliaSl^a' 2 

B lYM >BsS>,ee> ■ ■ ■ B ££ ',aa'/** (15) 

where 

-£>77',<5<5' = Cyjja 3 Cy'7)' ' a 3 Cr7<5a 4 Cr/' S' a t (1®) 

In fact, the matrix B is nothing but the ingredient of the 
denominator for po, "I"!*, and the right and left eigenvec- 
tor matrices arc introduced above as R and L and the 
eigenvalue matrix as p. We can write as B — RpL tr , and 
again using the property L tr ■ R = 1, we have, in the limit 
N -> oo, B*- 1 = Ro/^L^. We finally have, 

(AB) = {a2ai\AB\a l 1 a' 2 )-R, { 1 i)t (aa l )/p Q 

The numerical parameter used is, (3 = 10~ 6 . The case 
/3 = 10~ 7 gives negligible corrections to the entanglement 
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FIG. 3: The same as Fig. 2 for the local magnetic moment. 
The thick solid line (blue online) is the Bethe Ansatz result 
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FIG. 4: The ground state energy per site at half-filling as a 
function of the Coulomb energy U. From the top, the entan- 
glement n=l,2,3,5 (rectangles) and 7 (star, red online). The 
thick solid line (blue online) is the Bethe Ansatz result ||. 



n = 3 results below. The convergence criterion is \\C id ~ 
CeJI/IICJI < 5 • 10 ~ 5 - when this condition is met, 
the relative change in the largest eigenvalue [io often hits 
1CP 15 , the machine precision. 

Fig. 2 shows the ground state energy at U — 8 as 
a function of the electron concentration, n = 1 corre- 
sponding to half-filling. With the increase of the entan- 
glement n = 1,4,6 and 8, our result converges to the 
Bethe Ansatz result 17[. Fig. 3 shows the local mag- 
netic moment at U = 8 as a function of the electron 
concentration. Again, our calculation converges to the 
Bethe Ansatz result 17] . Fig. 4 shows the ground state 
energy at half-filling as a function of U . The results are 
from the top, n — 1, 2, 3, 5 and 7. There is a couple of % 
discrepancy at 17 < 1 from the Bethe Ansatz result 0. 
We have carried out the calculation for n — 9 and 10 for 
U = 0.01 (took about 10 hours using a single PC of about 
1 GHz processing speed) to get the ground state energy 
per site -1.25 and -1.255 to be compared with the ex- 
act one -1.2717. A rather slow convergence at U < 1 is a 
little surprise at first, but is understandable if we remem- 
ber that the kinetic energy term promotes electron itin- 
erancy, whereas the onsite Coulomb repulsion promotes 
electron localization. When purely itinerant, [7 = 0, the 
ground state is constructed by filling all the momentum 
states up to the Fermi level, giving the ground state en- 
ergy — 4/-7T at half-filling. If the free electron ground state 
is put in our form (|12[) . we would need a large entangle- 
ment number. On the other hand, when U — > oo, it is 
known from Bethe Ansatz that the ground state energy 
is at half-filling [17], which is just our result with entan- 



glement 1. In other word, the electron is fully localized 
and the mean-field treatment is good enough except its 
wrong, but not remotely wrong, prediction of antifcrro- 
magnetic long-range order, namely the Neel state. This 
is a delicate issue. In fact, in the limit U oo, the ID 
Hubbard model can be mapped onto the spin one-half 
antiferromagnetic Heisenberg chain which does not have 
long-range order. But the system is critical or quasi-long- 
range ordered in that its correlation functions fall off as 
a power of the distance 20]. In real materials, no truly 
ID quantum or 2D classical systems exist. There always 
exist 3D characters such as weak inter-chain or inter- 
layer couplings. Although weak compared to intra-chain 
or intra-layer interactions, these interactions are decisive 
for stabilizing the long-range order. 

In conclusion, the essence of the new method shall be 
summarized and possible future directions be discussed. 
First, the method is simple, general, not relying on ex- 
isting methods such as the cluster mean field theories 
and NRG. It only uses matrix algebra and fully imple- 
ments translational symmetry. Its application to other 
quantum systems in one dimension, namely quantum 
spins, bosons, and fermions with reasonable finite-range 
interactions and translational symmetry is immediate. 
Second, extension to thermodynamics with the use of 
standard procedure from quantum Monte Carlo, namely 
the quantum transfer matrix and its 90 degree rotation 
thereby reducing the thermodynamics to a similar eigen- 
value problem as treated in this paper, is straightforward. 
By switching between real and imaginary times, dynam- 
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ics should be handled as well. Third, and probably the 
most important and worth repeating the argument in 
the introduction, extension to the two dimension is also 
straightforward. In fact, mathematically, the extension 
from ID to 2D Hubbard models in our method should go 
similarly as in our study of the 3D Ising model based on 
the calculation of the 2D Ising model, the "Russian doll" 
structure. The only possible complication may arise from 
the anticommutation algebra in 2D fcrmions. At present, 
therefore, it would be safe to say that the extension to 
2D bosons and quantum spins is straightforward, but 2D 
fermions might need a further theoretical thought. 
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